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ABSTRACT 

The multifractal analysis of disorder induced localization-delocalization transitions is re- 
viewed. Scaling properties of this transition are generic for multi parameter coherent 
systems which show broadly distributed observables at criticality. The multifractal anal- 
ysis of local measures is extended to more general observables including scaling variables 
such as the conductance in the localization problem. The relation of multifractal dimen- 
sions to critical exponents such as the order parameter exponent /3 and the correlation 
length exponent V is investigated. We discuss a number of scaling relations between 
spectra of critical exponents, showing that all of the critical exponents necessary to 
characterize the critical phenomenon can be obtained within the generalized multifrac- 
tal analysis. Furthermore we show how bounds for the correlation length exponent V 
are obtained by the typical order parameter exponent Qfrj and make contact between the 
multifractal analysis and the finite size scaling approach in 2-d by relying on conformal 
mapping arguments. 

1. Introduction 

Anderson localization 1,2 has gathered interest for over three decades. Electrons 
in disordered conductors can undergo a transition to an insulating state. At zero 
temperatures the diffusion constant is a function of the Fermi energy ep and takes a 
finite value in the conducting phase while it vanishes in the insulating phase which 
is reached by crossing a certain value Ef = E c called the critical energy. This 
transition is due to a disorder induced localization-delocalization transition of the 
electrons wave functions. Localization occurs for strong enough disorder because 
of quantum interference effects brought about by the randomness of the disorder 
which is assumed to be static (quenched disorder). 2 

The modeling of disorder induced localization-delocalization transitions (LD 
transitions) refers to independent electrons of effective mass m moving in a ran- 
dom one-particle potential V(r), i.e. the Hamiltonian is H — 5— + V(r) where ir is 
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the kinetic momentum which may include magnetic fields. The potential is charac- 
terized by its mean value (V) and a finite range correlation function (V(f)V(f ')). 
A corresponding tight-binding version with on-site random energies and nearest- 
neighbor hopping is referred to as Anderson model. 

Transport properties of electrons are related to two-particle Green functions 
which are able to describe the transition probability from point r to point f ' at 
a given energy. 3 Here the notion of an n-particle Green function is introduced for 
the disorder average of a product of n Green functions involving retarded, G + , and 
advanced, G~ , Green functions. These are defined as the matrix elements of the 
resolvent operators corresponding to H, i.e. G ± {E) := (E — H ± ie)~ for e — > +0. 

For quasi 1-d systems (characterized by their length being much larger then their 
width) rigorous results are available 4 , saying that all states become localized as the 
length grows infinite. For localized states, the two-particle Green function shows 
exponential decay, \G + (E, r, f ')| cx exp — \r — f '\/£(E), from which a localization 
length £(E) can be identified. 

The lower critical dimension of the LD transition is believed to be d — 2. 5,6 While 
in the absence of strong magnetic fields or spin-orbit scattering most indications 5,7 
tell that all states are localized in the thermodynamic limit of 2-dimensional systems 
(though the localization length can be quite large), 2-d systems with strong magnetic 
fields undergo LD transitions which are believed to be responsible for the occurrence 
of the quantum Hall effect 8,9 (QHE) and are referred to as quantum Hall systems 
(QHS). There is now also striking numerical evidence that 2-d systems with strong 
spin-orbit scattering show an LD transition. 10,11,12,13 

In 3-dimcnsional disordered conductors the LD transition is believed to occur 
more generally 5 which has been confirmed numerically. 14,15,16,17 

Most studies on the LD transition (e.g. Refs. 14, 15, 17) were focused on the de- 
termination of the critical value of the Fermi energy, E Cl (or of the disorder strength) 
and on the critical exponent v of the localization length £. Since the late 80's, how- 
ever, it became clear that the wave functions at criticality are multifractal measures 
leading to a whole spectrum of critical exponents. 17,18,19,20,21,22 Nevertheless, only 
rare connection between these new exponents and the localization length exponent 
v was made. 17,20,23 

In this article, we wish to emphasize that the LD transition can be viewed as a 
prototype for a class of critical phenomena characterized by broadly distributed ob- 
servables at criticality. For this class the multifractal analysis provides a framework 
to investigate all of the critical exponents. Though, in the following, we often will 
refer to the LD transition for the sake of concreteness, we hope to make recognition 
of the generalities possible (a related approach has been put forward by Ludwig 24 
and by Ludwig and Duplanticr 25 ). 

The multifractal analysis is a scaling approach relying on the principle of "ab- 
sence of length scales" . Systems without characteristic intrinsic length scales (defin- 
ing the interrelations of the systems constituents) obey homogeneity laws with re- 
spect to rescaling. Let A be an operator or a complex valued function defined for 
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values x describing a length scale of descriptive nature (e.g. wavelength, system 
size) . Absence of length scales means that A shows a typical homogeneity law 

A(s ■ x) = s K ■ A{x) (1) 

where n is called the homogeneity exponent and s is a real number. In other words: 
the absence of length scales is reflected by the property that a rescaling of x can be 
compensated by a rescaling of the observable A. For real valued functions A(x) the 
solution of the homogeneity equation (|l|) is a power law: A(x) oc x K . Now, assume 
that A^ (x) is a functional of powers q of those observables which are involved in 
the definition of A{x), then the simplest situation that may appear is that n(q), 
defined by 

A^(s-x) = s Kiq} -A^ix), (2) 

is a linear function of q; the operators simply add their exponents. If n(q) shows 
a significant deviation from linearity we say that the scaling behavior of A{x) is 
anomalous. Multifractality will turn out to be a generic case of anomalous scaling 
behavior. 

The article is organized as follows. In Sec. 2 several aspects of critical phe- 
nomena, especially the notion of order parameters, correlation lengths and critical 
exponents, are reviewed and the LD transition is put into this perspective. It turns 
out that the LD transition shows unusual features, the most important of which 
is the occurrence of a dependence of the local susceptibility on the system size. 
In Sec. 3 we review the multifractal analysis of measures and the notion of f{a) 
spectra describing the scaling of broad distributions for local probabilities. Results 
for the wave function (defining a local measure) of QHS are discussed in detail. 
Generalizing to positive box-observables in Sec. 4 will allow us to relate the crit- 
ical exponent of the correlation length, v 1 to the multifractal exponents ao (the 
typical scaling exponent of box-probabilities) and X, the normalization exponent of 
a scaling variable. Furthermore, we discuss bounds for v in terms of multifractal 
exponents. In Sec. 5 we focus on correlations in local box-observables. This topic 
is touched already in previous sections, but is then treated more systematically. As 
a result we find scaling relations which determine all of the correlation exponents 
by the original /(a)-spectrum and the corresponding normalization exponents. By 
relying on conformal mapping arguments we conjecture scaling relations to deter- 
mine the multifractal spectra via finite size scaling methods in 2-d. Sec. 6 contains 
our conclusions. 

2. Aspects of Critical Phenomena 
2.1. Critical Exponents 

A defining feature of a critical phenomenon is the divergence of the correlation 
length £ c , defined by the spatial decay of the statistical correlations of a local field, 
{ip(0)ip(r)) oc exp (— r/£ c ), has to diverge at the critical point T c . In the sequel of this 
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article wc refer to the following nomenclature for critical phenomena. The physical 
states are described by state coordinates (we will refer to them by one symbol 
T) and by order parameters (m). For infinite system sizes the state coordinates 
exhibit a critical point T c where the order parameters vanish with a power law in 
t:= \T-T C \ 

m = (T < T c ) , m oc t (T > T c ) . (3) 

There exists a generating function F(t, h) where h is the conjugate field with respect 
to the order parameter m — — * 0) . In analogy to equilibrium systems we call 

F the free energy per unit volume. Higher derivatives of F also show power law 
behavior for h = 0, t — ► 0. For example 

d 2 F d 2 F 
X=Q^(h^0)*t-\ c =_(/ l ^0)oci- Q . (4) 

The exponents 7 and a (the latter may be negative) are called critical exponents of 
the susceptibility, \, and specific heat, c, respectively. 

To have a concept of spatial correlations we assume that the order parameter 
and its conjugate are local fields m(r), h(r) in a d-dimensinal system of volume V 
for which a generating functional 

Z (h(f), t) = J d[<p(f)] cxp {-H (Mr)},h(r);t)} (5) 

exists where TL is commonly called Hamiltonian although it docs not have to be 
hermitean. The functional F(h, t) is given as F(h, t) = — In Z(h, t) /V for the volume 
V — > 00. While the mean value m of the order parameter is given by the first 
moment of the generating functional, m = (<p(r)), correlation functions are given 
by higher order cumulants 
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(<p(r{)---<P(r-)) c :=- n7 ^ —(fc^Q). (6) 



which are related to moments in a combinatoric way, for example: (f{r)) c = {fir)) 
and 

(V(r>(f '))c = MrMr')) - W')) =■ x(\r~ ?'\) ■ (7) 

x(r) is called local susceptibility, since the generating formalism allows to conclude 
/ d d rx{r) = X- The correlation length £ c can be identified from x(r) cx cxp— r/£ c 
and is assumed to diverge with a power law approaching t = 0, 

&(t) cx t~ v , (8) 

where v is called critical exponent of the correlation length. 

To illustrate the analogies of the LD transition with critical phenomena as out- 
lined here we refer to the fact 3 that the LD transition is characterized by the analytic 
behavior of the disorder averaged two-particle Green function, 



9 (\r~r 'I) = < r\G+\r> >< r'\G~\r>, (9) 
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with respect to the infinitesimal parameter e in G ± = (E — H ± ie)^ 1 which distin- 
guishes retarded and advanced Green functions and controlls the long time averages 
of transport quantities such as the diffusion constant D(E), 

p(E)D(E) = lim -\e f d d rr 2 g 2 {r), (10) 

where p(E) = ^< r\G~(E) - G+(E)\r> is the average density of states. Thus, 
it seems natural that a critical phenomenon description of the LD transition may 
start from a generating function for disorder averaged Green functions. For our 
purposes we do not need to go into technical details and leave the explicit form of 
the corresponding Hamiltonian open. We only like to refer to some of the main 
features of such theories (for details see e.g. Refs. 26, 27) which can be summarized 
as follows. 

The degrees of freedom are matrix fields Q with a block structure Q pp where 
p,p' G {±} correspond to retarded (+) and advanced (— ) Green functions, respec- 
tively. The source field is the infinitesimal parameter e and the corresponding source 
term in the Hamiltonian is e (Q ++ — Q~~) leading to 

5 -^- = (Q— ) - (Q++) = < r\G-\r >-< r\G+\r > (11) 

which is proportional to the average density of states p(E). Thus, the density of 
states in disordered electron systems seems to be a candidate for an order parameter 
of the LD transition. Unfortunately, the density of states is a smooth function of 
energy 3 ' 28 unable to reflect the LD transition which is characterized by the vanish- 
ing of the diffusion constant related to the susceptibility via Eqs. [h]). Before 
coming back to this peculiarity of the LD transition let us proceed in listing common 
knowledge about scaling ideas applied to critical phenomena. 

2.2. Scaling Relations 

The crucial assumption of any scaling approach to critical phenomena (see e.g. 
Ref. 29) is that the critical exponents a, (3, 7 and v have their origin in the diver- 
gence of one relevant length scale, the correlation length £ c . Since thermodynamic 
quantities are generated by the free energy we make the homogeneity assumption 
that the free energy, in the vicinity of the critical point, can be written as 

F(t;h) = at a f ± (h/t A )+f r (12) 

where f± is a regular function of the argument h/t A and f r stands for regular terms 
which do not show scaling behavior for h — * 0, t — > 0. By differentiating F, the 
following scaling relations 

2/3 + 7 = 2-a, a = 2- a (13) 

can be concluded. The scaling exponent v of the correlation length £ c which relates 
the divergence of £ c with the vanishing of t is related to these exponents by the 
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following scaling argument. Consider a finite system £1 of volume L d for which the 
free energy per unit volume is 

F&h) = (14) 

A change of L in Z^ltojho) to L' = XL for fixed values to, ho is expected to be 
compensated by an appropriate change in t, h such that the partition sum remains 
unchanged, i.e. 

Z L (t ,ho) = Z L ,(t',h'). (15) 

This procedure defines implicitly functions t' (L, L'), hi (L, L') with t(L,L) = to, 
h(L,L) = ho- Consequently, 

F L (t ,ho)=F L >(t',ti)-\ d . (16) 

Choosing the value L — £ c (to) we see F^u^ito, ho) — FL>(t r , h') ■ £r d (io) and find 
a scaling relation 

U = d-v. (17) 



We can cast the scaling relations Eq. (|13|) in a form, where only a, (3, 7 and v 
appear 

2(3 + 1 = 2-a = d-v. (18) 

This means that only two of these four exponents are independent. 

We introduced the correlation length with the help of the local susceptibility. 
At the critical point, where £ c diverges, the correlation function can show a homo- 
geneity law * 

X(r) oc t~~ z . (19) 

Owing to the general relation between the global susceptibility x an d the correlation 
function x( r )i X — J d d rx{r), and by introducing the correlation length £ c (i) as a 
cut-off length in the integral over x(r), which diverges at t = 0, we find another 
scaling relation which is independent of the homogeneity assumptions for the free 
energy 

7 = v (d-z) . (20) 

Together with the scaling relations of Eq. ( jl8| ) this tells that only two of the five 
exponents a, [3, 7, v and z are independent. 

We now come to a crucial point for what follows on the multifractal analysis 
of broadly distributed observables at criticality. Imagine that in addition to the 
distance r the finite system size L is in a regime of absence of length scales. Assume 
that the function C(r, L) := (ip(r)(p(f ')) L shows a homogeneity law with respect to 
both lengths, r and L, 

C(r, L) cc r~ z L~ v . (21) 

*Note, that in the literature on critical phenomena the exponent z is often written in a form 
z = d — 2 + r\. We will avoid this notation since it could lead to confusions in the context of the 
LD transition where an exponent r) is used sometimes in a different context (see the discussion 
below). 
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The exponent y describes the system length dependence and z ^ z describes the 
distance dependence. That such a situation will indeed appear in the LD transition 
problem will be discussed in Sec. 2.4 and in more detail in Sees. 3.3, 5.2. A similar 
reasoning now shows that 

y=v(d-(y + z)) , y + z = z (22) 

where the last relation follows from consistency with Eq. (pp|). Notice that identi- 
fying the exponent n in z — d — 2 + 77 as the correlation exponent means that such 
rj is not determined by 7 and v alone, but requires the additional knowledge of y. 

2.3. (3-Functions 

A method to test scaling assumptions and calculate critical exponents is the so-called 
renormalization group. The definition of the renormalization group is not unique, 
but can in general be described as a transformation acting on the Hamiltonian Ti. 
of the system where a length scale serves as transformation parameter. Instead 
of introducing a finite system size, one can also work in the thermodynamic limit 
considering the system to be defined on a lattice. Changing the lattice constant 
clq to Lb = bao, and changing the coordinates t 1 h in an appropriate way, one tries 
to keep In Z fixed. After many iterations of this procedure (which defines the 
renormalization group) one ends up with a new Hamiltonian the structure of which 
should not have changed except for the change in the coordinates and that the fields 
are now defined on a new lattice. That such a procedure indeed works makes the 
model renormalizable. Otherwise one had to introduce a number (in the worst case 
an infinite number) of new terms and coordinates into the new Hamiltonian in order 
to keep \nZ fixed. Thus, the application of the renormalization group to the free 
energy density leads to 

F := F[t, h) = b- d F(t', ti) (23) 

where t' = b Vt t 1 h! = b Vh h in the critical regime and yh, yt are scaling exponents. 
Choosing the factor b such that t' = const., we arrive again at the scaling relations 
of Eqs. (|l3],[l8]). By considering the system size L being the scaling parameter, 
an extension of renormalization group ideas can be made, based on the following 
assumption. Far away from the critical regime (where power law scaling holds) there 
still exist scaling variables A(L) which fulfill a restrictive functional equation: 

A(bL) = f (A(L),b) . (24) 

For simplicity we consider the case of only one scaling variable here. This means 
that the value of A for system size bL only depends on the value of A for system 
size L and the scaling factor b where / is called a scaling function. This assumption 
is very strong since, in general, A will depend on many microscopic parameters. 
The idea behind this assumption is that most of these parameters will become 
unimportant for describing the flow of the scaling variable with increasing system 
size, and that the scaling function / is the only information needed to determine 
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this flow. The existence of a scaling function / is equivalent to the existence of a 
/3-function 

which is indeed a function of In A alone. Provided the /3-function is smooth, the 
flow is given by the solution of the differential equation 

„ ., . . din A . „. 

/3 0bA) = — . (26) 

Eq. (^6|) is called a renormalization group equation. The regime where /3 can be 
linearized around a fixed point A* of the flow {(3 (In A*) = 0) is called critical regime, 

/3(lnA) =0 • (InA-lnA*) . (27) 

Here (3' is the slope of the (3- function at A*. Starting from a system of size L = Lq, 
with A = Ao chosen close to A*, and turning on the renormalization flow until the 
system reaches a size £ c marking the bound of validity of Eq. (27) (which is here by 
definition the correlation length), one finds 



/InA-lnA* 



L J VlnAo-lnA* 



(28) 



The critical regime is narrow, we expand A(£ c ) and Ao around A* and get for the 
correlation length 

/ A - A* \ 1//3 ' , N 

& = Lai - A ■ (29) 



> -A*. 

Consider the situation where the width of the critical regime AA = Ao — A* is 
triggered by the parameter t (leaving Lq and A — A* unchanged), the divergence 

£ c oc (AA)"" oc (30) 

follows, where the exponent v is given by the inverse slope of the (3- function at the 
fixed point, i.e. 

V =l/I3'{k*). (31) 



2.4. Unusual Features of the LD Transition 

It is now time to put the LD transition into the perspective of the critical phenomena 
terminology. A common feature of field theoretic approaches to the LD transition is 
that disorder averaged n-particle Green functions are generated by a field theoretic 
generating functional in the sense of Sec. 2.1. As already mentioned, the density of 
states which is a one-particle Green function appears as the formal order parameter 
and the conjugate field to the density of states is the infinitesimal parameter of 
the Green function. The density of states does not show the LD transition, but 
is a smooth function of the energy which, in the problem of the LD transition, is 



Multi-racial Analysis ... 9 



the analog of the temperature in equilibrium phase transitions. Consequently, the 
"critical exponent" of this formal order parameter is zero 

= 0. (32) 



The LD transition cannot be described by the formal order parameter, Eq. (|l] 
of the generating functional. The LD transition can only be obtained by analyzing 
quantities which are related to the two-particle Green function such as the diffusion 
constant D(E). 

Furthermore, the scaling relations, Eqs. ( |l8| , po|) , tell that the correlation expo- 
nent of the local susceptibility (a two-particle Green function) vanishes, i.e. 

5 = 0. (33) 

This result seems to rule out a power law behavior of the two-particle Green function 
(e.g. the density correlator). However, Chalker and Danicll 30 demonstrated, by 
numerical calculations, the density correlator to show a power law behavior at the 
LD transition of a quantum Hall system, oc 7- -0 - 38 . I n addition, by relying on an 
exact inequality obtained by Chalker 31 one has analytical evidence that the density 
correlator at the critical point cannot be described by a vanishing distance exponent. 
This phenomenon is often referred to as anomalous diffusion. Fortunately, from 
our discussion of correlation exponents in Sec. 2.2 where we distinguished between 
different types of distance exponents z and z we can fairly conclude that the value 
of 0.38 does not correspond to z, but corresponds to the distance exponent z in a 
regime where the correlation function shows power law also with respect to another 
length scale which can be interpreted as a system size L being much smaller than 
the correlation length £ c , i.e. 

z = 0.38^ 5 = 0. (34) 

The phenomenon of anomalous diffusion and the unability of the density of states to 
describe the LD transition are two of the main unusual features of the LD transition. 
By Eq. (|34|) it is indicated that both are not independent of each other. 

Another unusual feature of the LD transition which, at first sight, seems to 
be independent of the previous ones was observed by Altshuler et al. , 32 ' 33 They 
reexamined the phenomenological scaling theory for the LD transition of Abrahams 
et al. . 33 In this scaling theory the conductance g{L) of a cube with linear dimension 
L is considered to be a scaling variable in the sense of Sec. 2.3, i.e. g obeys a 
differential equation in terms of a /3-function 

It! -'<-"> < 35 » 

which is a unique function of In g. As shown by Altshuler et al. 32,33 the conductance 
of mesoscopic systems highly depends on the individual properties of a given system 
(e.g. on the concrete disorder potential) . As a result an ensemble of different systems 
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reflects a broad distribution in g which cannot be characterized by the mean value 
of g alone. Consequently, the mean value of the conductance {g) is not a suitable 
scaling variable. Instead, in any scaling approach to the LD transition one has to 
consider the whole distribution function of the conductance. However, it may be 
possible to apply the ideas of the renormalization group to certain parameters of 
this distribution function. Following Shapiro 3 we refer to an s-parameter scaling 
theory for the LD transition if s parameters <7 re ii ... s obeying scaling equations 

din g ro i„ 

-^=^{{9^k}k=i,...,s) (36) 

are required to determine the distribution function of a scaling variable g. If, at least 
in the vicinity of the transition point ({g* elk }k=i,...,s — 0)), only one relevant 
length scale £ c and one type of scaling variable g Te i{L) exists, then the critical 
exponent v of the correlation length £ c is defined by 

g iel (L) - 5r * el oc for L < £ c (37) 

and l/v is given by the slope of the corresponding (3- function at g* cl . A context 
where this concept will serve to be fruitful is provided by the generalized multifractal 
analysis to be discussed in Sec. 4. There we will also see that anomalous diffusion 
and broad distributions of physical quantities are deeply connected. 

3. Multifractal Analysis of Measures 

After having discussed some aspects of critical phenomena, where the absence of 
length scales is reflected by the existence of critical exponents, the multifractal 
analysis of measures is reviewed in this section. 35 ' 36 ' 37 ' 38 This analysis is appropriate 
to describe self similar local observables which can be interpreted as measures. It 
can also be viewed as an extension of the fractal dimension approach to self similar 
geometric objects invented by Mandelbrot. 39 

3.1. Scaling of Moments 

To introduce the notion of the multifractal analysis and for the sake of concreteness 
we consider a quantum Hall system at criticality. For a finite two dimensional 
quantum Hall system of linear size L with double periodic boundary conditions we 
study the electrons wave function ?p(r), the modulus of which defines a normalized 
measure. The probability for an electron to be found in a box of linear size Lb is 
given by the box-probability 

P(L b ) := J d 2 r\m\ 2 - (38) 

box 

Covering the system by a mesh of N(Lb,L) boxes the fractal dimension D of the 
wave functions support is defined by N(\) oc \~ D where N(X) is the number of 
boxes with non- vanishing box-probability for a given ratio X — Lb/L. Since the wave 
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function is never exactly zero in any box, N(X) equals the total number of boxes 
N(Lb, L) and, consequently, the fractal dimension is the Euclidean dimension of the 
system, i.e. D = d = 2. At the level of the fractal dimension D the wave function 
shows no interesting fractal behavior. Let us now focus our attention to the scaling 
behavior of the box-probability. The normalization condition X)f^ Pi(Lb) = 1 
yields the scaling behavior for the average 

(P(L b )) L *\ D . (39) 

Here the average of quantities Ai (corresponding to the box numbers i) is defined 
by {A) L :— w^y Si A* ■ Let us consider higher moments of the box-probability 
P{Lb) with respect to this averaging procedure. The general assumption underlying 
the multifractal analysis is that for a finite interval of values of A the moments 
{P q (Lb)) L show power law behavior indicating the absence of length scales in the 
system, i.e. 

(P q {L b )) L cx X D+T ^ (40) 

where r(q) does not depend on A. In quantum Hall systems the assumption about 
the absence of length scales means that on the one hand the correlation length £ c 
of the LD transition has to be much larger than the lengths L b ,L. On the other 
hand microscopic length scales Z, such as the cyclotron radius r c , have to be much 
smaller than L b , L. In summary the condition for the multifractal analysis to be 
useful is 

I < L b < L < £ c . (41) 

For finite systems the states which are candidates for the multifractal analysis are 
states with localization lengths £ being much larger than the system size. This 
corresponds to the critical states of the LD transition in finite systems. In the 
thermodynamic limit L — > oo such states can only be found at the critical energy 
E c . In this limit A goes to zero and the function r{q) can be defined uniquely for 
the critical state by 

T(9):= lTo hTA D - (42) 

For finite systems one can give an estimate for r(q) by considering the slope in a 
plot of In ((P q (L b )) L ) versus In A over an interval of values A which fulfill the above 
mentioned constraints. Now we want to discuss how multifractality is reflected by 
the properties of the function r(q). Note that by the construction of (P q (L b )) L wc 
have found an observable in the sense of Eq. (||) . The function r(q) corresponds to 
the homogeneity exponents n(q) and multifractality of the wave function means that 
r(q) is a non-linear function of q. According to the normalization condition, t(1) 
, t(0) = —D and r(q) can thus be parametrized by the generalized dimensions 
D(q), defined by 

r(q) =: D{q) ■ (q - 1) , (43) 

and the fractal dimension D — D(0). Consequently, a single- fractal wave function 
is characterized completely by the fractal dimension D, i.e. D(q) = D. D(q) being 
not a constant function distinguishes multifractals from single-fractals. 
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The reason why the multifractal analysis is a powerful method to analyse scaling 
behavior comes from the fact that one can make very general statements about the 
analytic behavior of the r(q) function and the D(q) function, respectively. Assuming 
smoothness of both functions for q being real numbers one can derive the following 
results: 38 ' 40 

• The function r(q) is a monotonically increasing function with negative curva- 
ture, i.e. 



^>0,^<0. 



dq 



' dq 2 



(44) 



The generalized dimensions D(q) are positive, monotonically decreasing and 
bounded by D± oa := D(q — > ±00), i.e. 



dD(q) 
dq 



< , < L>oo < D(q) < D. 



(45) 



• The function r(q) has asymptotically constant slopes given by D± oa 

Here a sketch of the proves of these statements is given (for details see Refs. 38, 40) 
The main part of the derivations relies on the homogeneity, reflected by power laws, 
already for finite ratio A, and on the normalization condition which keeps the box 
probabilities less or equal to unity. First, the monotonicity of r(q) is a simple 
consequence of Pi < 1 and, thus, J2i^i < -^f for q > q' . A much stronger 
statement is the monotonicity of D{q) which can be shown by writing D(q) as 



9-1 



1/(9-1)' 



(lnA)" 



(46) 



with the help of a generalized Holder inequality, 



l/r 



> 



l/r' 



for r > r' . 



(47) 



Notice that for distinct values of q and q' the equality D(q) = D(q') only holds in 
the case of a constant function D(q) = D. For each ratio A there exist maximum 
and minimum values for the box-probabilities. Now raising the box-probabilities to 
positive powers q — > 00 means that in the sum ^ P? only the maximum values 
contribute significantly yielding the dimension > by 



c (L b ) cx A 1 



(48) 



A similar consideration for q — > — 00 shows that the minimum values of the box- 
probabilities dominate in the sum yielding the dimension D_oo by 



P min (L fc ) cx X D - 



(49) 
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Thus, the r(q) function has asymptotic slopes given by D± oa which shows that 
the limits of D(q) for q — > ±00 exist and are just given by D± oc . Notice that we 
referred to finite ratios A guaranteeing the existence of maximum and minimum 
values. Extrapolating from small but finite values of A to the thermodynamic limit 
A — > becomes possible due to scaling, i.e. the absence of a length scale in the 
system. To complete the proof we have to show that r(q) has negative curvature. 
Due to the smoothness assumption 

exist and are continous functions of q. The definition of r(g) yields 

- P'ifr* -° 2 '-'))-'° A 

. ^.SiiS^^aXo (51) 

Z(q, A) 



where 



Z(g,A):=^P? (52) 



is called the partition sum of the multifractal. The equality in the relation a' (q) < 
only holds iff a(q) = D, i.e. in the single-fractal situation. From Eq. ( |50| ) one can 
see that the objects characterizing multifractals are constructed in a similar way 
to thermodynamic quantities where the partition sum Z(q, A) is replaced by the 
canonical partition sum Z((3, V) in Boltzmann statistics. Indeed, one can write 
down for each of the quantities like q, lnA,r(q) ... its thermodynamic counterpart 
(see e.g. Ref. 41). Keeping this analogy between multifractals and thermodynamics 
in mind one can imagine that a violation of the smoothness condition for r(q) 
could be viewed as an intrinsic phase transition in each multifractal wave function. 
Such kind of phase transition (which has not been observed in systems showing LD 
transitions up to now) should not be confused with the LD transition itself. 



3.2. Distributions 

The r{q) function describes the scaling behavior of moments, i.e. of (P q ). We 
wish to describe the whole distribution function U(P; A) that corresponds to these 
moments. For a function F(P; A) the corresponding average value reads 

1 

(F(P;X)) L = J dPTl(P;\)F(P). (53) 
o 
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Since we are interested in the scaling behavior of the distribution function itself we 
change the variable P to a by defining 

lnP 

a := in! 

fl(a;\)da := U(P;\)dP. (54) 
Now, average values can be calculated with the help of II(a; A) via 

oo 

(F(P; X)) L = J da fi(a; X)F(X a ) . (55) 
o 

Absence of length scales forces the distribution function II to show power-law scaling 
with respect to A. Since, in the thermodynamic limit A — » 0, the values r(q) are 
uniquely defined one can apply the method of steepest decent when calculating the 
average for F = P q . This leads to the following conclusion 

• There exists a positive function, the /(a)-spectrum of the multifractal, char- 
acterizing the scaling behavior of the distribution function IT, and is related 
to r(q) by a Legendre transformation. 

fi(a;A) oc \-f( a )+ D 

f(a(q)) = a(q)-q-T(q). (56) 
Here the function a(q) is the same as introduced in Eq. (j50|). 

Using the analogy between multifractality and thermodynamics one can check that 
f(a) corresponds to the microcanonical statistics and resembles the entropy as a 
function of energy. Let's focus on the meaning of the function f(a) in the context of 
multifractality. This function describes the scaling behavior of the whole distribution 
function of box-probabilities in the absence of length scales. In the single-fractal case 
the function f(a) shrinks to the point (P/, D) in an (a, f(a)) diagram. Consequently, 
the distribution function II is singular, and the distribution function with respect to 
the variable P, II(P; A), is a narrow distribution on all length scales. In the single- 
fractal case II(P; A) is typically a gaussian distribution and can be characterized by 
a few cumulants. 

In the multifractal case there appears a new and unexpected behavior of the 
distribution function II(P; A) which we want to describe in the following. Due to 
general properties of Legendre transformations we have 

• f(a) is a positive, single- humped function of negative curvature on a finite 
intervall of a values: 

D 00 <a<D_ 00 , 0<f{a)<D, f(a )=D (57) 



where uq = a(q = 0) > D. 
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• /(a) terminates at the points (-D±oo,0) with infinite slopes and has slope 
f'(a) = 1 at the point (a(l),a(l)) where a, as well as f(a), equal D(l). 

These statements can most easily derived by writing the functions r(q), a(q) 
and f(q), defined by f(a(q))=: f(q), as 42 

r(q) = ±-\nZ(q,\) 
a{q) = ^J^Vf^Pi 

i 

/(<?) = j^$Z^' ln ^- ( 58 ) 

i 

Here the g-dependent normalized quantity 

^(q,L b ) := J , ^ b) (59) 

is a generalization of the original box-probability Pi and is called the q-microscope 
since it increases the large box-probabilities for positive values of q and increases the 
small box-probabilities for negative values of q, respectively. The second of these 



equations ( 58 ) is only a reparametrization of Eq. ( |50| ) and by Legendre transforming 
r(q) one can easily check that f(q) is correctly described by the third equation of 
(|58|). We mention that numerical calculations of /(a)-spectra are of higher precision 
when using Eqs. ( |58| ) compared to numerical Legendre transformation of r(q). 

Let us summarize the instruments for describing the multifractal scaling behav- 
ior of normalized box-probabilities: the r(q) function describes the scaling behavior 
of moments and f{a) describes the scaling behavior of the corresponding distribu- 
tion function. Both of them are related by Legendre transformation. The functions 
D(q) and a{q) can also serve to describe the multifractal nature of box-probabilities 
but have, for our purposes, less appealing interpretations. 

Motivated by the single-humped shape of the /(a)-spectrum we try to give a 
physical interpretation of multifractality, at least in the context of the LD transi- 
tion. It is illuminating to approximate the /(a)-spectrum by a parabola and to 
see to which kind of distribution function the parabolic approximation (PA) corre- 
sponds. The PA, which fulfills most of the desired constraints, is determined by one 
parameter, ao, besides the geometric fractal dimension D, 20,37 i.e. 

This corresponds to a log-normal distribution of the box-probabilities 

* (a;A) ^l^- lnA }- (61) 
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A log-normal distribution is the prototype of a broad distribution which cannot be 
characterized by a few cumulants. For example, the n-th cumulant grows exponen- 

2 

tially with n, oc e" . Scaling requires that the distribution is broad on all length 
scales and with the normalization condition on the box-probabilities the form given 
by Eq. ( |sT[ ) can be viewed as the paradigm of a distribution function which is broad 
on all length scales. The physical interpretation of such a behavior is the following. 

In phase coherent disordered conductors the actual value of a local box-quantity, 
like Pi(Lb), depends on a large number of conditions simultaneuosly. This can be 
simulated by writing Pi(Lb) as a product of a large number of independent random 
factors, 

Pi(L b ) =p - pi -p2 ■ ••■ , (62) 

and the central-limit-theorem tells that one can expect Pj(Z&) to be log- normal dis- 
tributed. The scaling condition (absence of length scales) requires the distribution 
function to exhibit power law behavior. Thus, we conclude 

n(^ ;A)ocA S-/(lnP/lnA) (63) 

in A 

with /(a) being almost parabolic. For completeness we list the PA for D(q) and 
a(q) 

D(q) = D-q{a Q -D) , a(q) = a - 2q (a - D) . (64) 

However, the parabolic approximation can never be exact since it violates some 
of the constraints we have derived for /(a); mainly positivity and the boundary 
conditions at (D± OQ , 0). PA can only serve as a good approximation in the vicinity 
of the most probable value ao for the scaling exponents a. The breakdown of the 
PA can be indicated by either the values a ± = ao ± 2yD[otQ — D) (where f{a) 
vanishes in the PA), or by q + = 2 (c^-D) ( wnere T~(q) looses monotonicity in the 
PA). a± give rough estimates of D too . 

The most probable value of the box-probability corresponding to ao is given by 
the typical value of P(Lb), 

P typ (L b ) := exp{(lnP(i 6 )) L } c< X a ° . (65) 

which is a geometric type of mean value. 

Finally, we present results of the multifractal analysis for wave functions in the 
critical regime of a finite quantum Hall system. 20 In Fig. |l| the squared amplitude 
of a wave function is shown where increasing darkness correponds to higher prob- 
ability to find an electron at the corresponding point in the system. The picture 
demonstrates the self-similar structure of the wave functions and the tremendeous 
amount of fluctuation of local probabilities, even within one given state. In Fig. || 
f(a) spectrum for this wave function is graphically depicted, showing that the box 
probabilities are almost log-normal distributed and show power law behavior. The 
most interesting values of scaling exponents are 

a = 2.3 ± 0.07 , D x = 0.95 ± 0.1 , = 3.7 ± 0.1 . (66) 
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Fig. 1. Squared amplitude of a multifractal wave function in the critical regime of a quantum 
Hall system. Increasing darkness corrcponds to higher probability of finding an electron (on a 
logarithmic scale). 



Fig. 2. The f(ce) spectrum of a multifractal wave function at the center of the lowest Landau 
level. 
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Fig. 3. Comparison of /(a)-spectra for three different wave functions in the critical regime of a 
quantum Hall system. The dots correspond to the /(o;)-spectrum and other curves correspond 
to the parabolic approximations for the /(a)-spectra of this and two distinct states. 

The fractal dimension is, of course, D = 2. 

Striking is the observation that the f(a) function seems to be universal for all of 
the critical states. This can be seen in Fig. || where the /(a)-spectrum of one state 
is shown together with the parabolic approximations of this and two distinct states 
choosen in the critical regime. The coincidence is remarkable. Calculations of 
/(a)-spectra, carried out by Pook 20 for the same model system, but with different 
realizations of the random potential and different system sizes, confirm this result. 
Further evidence was provided by calculations of Huckestein et al. 22,43 for a tight- 
binding version of the quantum Hall system which give the same values for the 
quantities ao, £*±oo within the error bars. This means that the LD transition of 
quantum Hall systems (at least for the lowest Landau level where all of the numerics 
has been performed) is uniquely characterized by those critical numbers. 

3.3. Inverse Participation Numbers 

In this subsection we describe how disorder-averaged Green functions are involved 
in the multifractal analysis via inverse participation numbers. Thereby two appli- 
cations will be adressed: the first concerns the calculation of the critical exponent 
v of the correlation (localization) length in a quantum Hall sysytem. The sec- 
ond concerns the anomalous diffusion found by Chalker et al. 30,44 as mentioned in 
Sec. 2.4. 

The role of fractality for the LD transition was already observed by Aoki 45 in 
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1983 by looking at the inverse participation number defined by 

V = j d d r\iP{f)\ A (67) 
n 

where f2 denotes a (i-dimensional system with linear dimension L. He gave an 
argument that the wave function tp(r) of electrons at the LD transition point shows 
fractal structure. If the wave function is uniformly distributed, V oc L~ d , the 
participation ratio p = {VL d )~ x is constant. In the localized regime V w and 
p vanishes in the thermodynamic limit. At the transition point where the wave 
function is extended the participation ratio still vanishes in the thermodynamic 
limit. Consequently, V scales with a power d* < d. By a box-counting method 
Aoki 46 calculated the fractal dimension of wave functions of a QHS at the transition 
point to be approximately 1.5. Although his method was not free of systematic 
errors, plots of the electron density |^(r)| 2 indicated self-similarity. 

Castellani et al. 19 considered the a d = 2 + e" expansion of Wegncr's non linear 
cr-model 47 for the generalized inverse participation numbers 

V k\ = J d d r\iP(r)\ 2q oc IT 7 "* <«> (68) 
n 

and concluded that the wave functions at the transition point show multifractal 
behavior since r*(q) is not equal to d(q — 1) but defines a set of fractal dimensions 
d*(q): 

r*(q)=d*(q)(q-l) . (69) 

To show the relation between r*(q) and r(q) we consider "P^l calculated in a tight 
binding model with lattice constant b, 

V^=J2mi)\ 2q = Z(q,X = b/L). (70) 

i 

If both quantities b and L are in a regime of self-similarity the identity of r* (q) 
and r(q) follows immediately. However, they are different in the continuum limit 
(b — » 0) for fixed L since self-similarity breaks down once b is smaller than the 
microscopic length scales I. On the other hand, in the continuum limit, we can 
consider Z(q,X = b/L) to be an approximation for which becomes more and 
more accurate as L goes to infinity. If self-similarity is preserved up to length scales 
L — > oo, T*(q) — > r(q) can be concluded. The assumption of self-similarity for 
L — > oo is fulfilled at the critical point of the LD transition where the only relevant 
length scale, the localization length £, diverges. Thus, at the critical point of the 
LD transition 

T*(q)=T(q). (71) 

On the other hand the correlation length £ c is defined as an upper bound for power 
law behavior and d*(q) can be determined from the scaling behavior of with 

P w «^ r * (?) . (72) 
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The critical exponent v of the correlation length £ c is defined by 

i c oc A"" (73) 

where A is the difference of a critical observable (e.g. the energy) from it's critical 
value. In analogy to this, Wegner 47 introduced critical exponents 7r(g) for the 
generalized inverse participation numbers 

pWl „. . (74) 

Thus, we conclude a scaling relation between multifractal scaling numbers r(q) and 
critical exponents ir(q) and i/ 17,48 , 

n(q) = v ■ r(q) , (75) 

and Tr(q) depends on q in a non-trivial way if the wave functions are multifractals. 
Still, Eq. @ doesn't yield v as a function of only multifractal scaling numbers. 
Nevertheless, we can see what was overlooked by Hikami 49 when calculating the 
critical exponent v for a QHS (restricted to one Landau level) to be 1.9 ±0.2. What 
he did calculate was the exponent 7r(2). He found 7r(2) = 3.8 ± 0.4 The result for v 
given above follows if P' 2 ' is assumed to scale with the localization length £ as £~ 2 , 
ignoring that multifractality leads to the anomalous scaling behavior oc £~ D ( 2 \ 
Taking multifractality into account, D(2) = 1.62 ± 0.02, 50 we get 

i/ = 2.4±0.3. (76) 

This is in agreement with the high accuracy result, v = 2.34 ± 0.04, of Huckestein 
and Kramer 51 obtained by a finite size scaling technique. 

As a further application to the QHE we show that multifractality relates the 
exponent r\ of anomalous diffusion to D{2). We have already anticipated that the 
spectrum of multifractal dimensions has universal features for states in the vicinity 
of the LD transition since we omitted to distinguish between inverse participation 
numbers of individual states and their ensemble average. In fact, as already dis- 
cussed, different wave functions of several systems show, in the critical regime, the 
same /(a)-spectra within the error bars. 

The observation of anomalous diffusion (see Sec. 2.4) means that the averaged 
two-particle Green function 



\G+(r,E)\ 2 = |(f|(£-fr + ie)-i|0)| 2 (77) 

behaves as 



\G+(r,E)\ 2 cxr-" +2 - d (78) 

This power law, with a value rj = 0.38 ± 0.04, was found 30 in quantum Hall systems 
(d=2) for length scales r in the regime of multifractality, i. e. / -C r <C £• Replace 
the ensemble average of the inverse participation ratio ( |69| ) of wave functions ip Sa (i) 
with 

(25) := J2 E 1^- Wl 29 • S ( E ~ e «VE § ( E - £ «) ( 79 ) 
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where i^ Sol (i) corresponds to energy eigenvalue e a . The ensemble averaged inverse 
participation number V(E) is now defined with respect to energy. It can be calcu- 
lated from the Green function 47 



where p(E) = lim - J d d r |G+(r, E)\ 2 is the density of states. The expected scale 

independence of this quantity is excellently confirmed by numerical results 20 . Thus, 
with the ansatz 



one concludes 



and consequently 



\G+(r,E)\ 2 =1 {E,L)r~ z (81) 
V(E) oc j(E, L) oc L z ~ d , (82) 
r 1 + d-2 = z = d- D(2) . (83) 



More precisely, in the case of multifractal behavior one expects the Green function 
to have a spectrum of exponents and the one occurring in Eq. ( |S3|) should be inter- 
preted as the dominating one. In the regime of mult ifr act ality, besides the distance 
r, a second length scale, the system size L, appears in the correlation function \G + \ 2 
(Eq. (|8l|)). This was anticipated in Sec. 2.4 when discussing the difference between 
the distance exponents z and z. Using the terminology of Sec. 2.2 we identify the 
exponents y and z to be 

y = -z = D(2)-d , 5 = y + z = 0. (84) 

In a single-fractal situation, z = and y = z. Thus, in systems without multifractal 
structure there is no need for introducing two exponents when discussing correlation 
functions in finite systems. The numerical value 50 D(2) = d*(2) — 1.62 ± 0.02 is in 
accordance with our considerations, and anomalous diffusion can be interpreted as 
being a direct consequence of the multifractal nature of wave functions at the critical 
point of the LD transition. This was already conjectured in Ref. 30. There is no 
contradiction between a finite value of 77, as interpreted via the distance exponent 
z instead of z, and the vanishing of the exponent /3 corresponding to the formal 
order parameter of the LD transition. The physical reason for the phenomena 
of anomalous diffusion and of mult ifr act ality is the quantum coherence of highly 
disordered systems. 

4. Genralized Multifractal Analysis 

So far, only measures, i.e. box-probabilities P, have been considered. To get into 
contact with usual critical exponents we introduce a generalized multifractal anal- 
ysis for non-normalized observables. 20 A related approach in the context of field 
theory was put forward by A. Ludwig. 24 
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Consider a physical observable Q which can be defined for systems of any linear 
extension L. For simplicity, we assume the system to be a box of volume (Lb) d 
where d is the euklidean dimension of the system. Thus, for any of such boxes 
we have a box-observable Qi(Lb) similar to the box-probability in the foregoing 
sections. Examples are the magnetic moment or the two-terminal conductance of a 
rf-dimensional cube with volume (Lb) . 

We are going to study the statistical and scaling properties of box-observables. 
We restrict to positive values for Q, 

Q(L b )>0, (85) 

to avoid cancellations in averaging procedures and to apply the results of the fore- 
going sections. At first glance, this seems to be a strong restriction. However, 
still a large class of observables can be studied which is specified by the following 
consideration. Let a physical observable A(Lb), measured on a box, be an element 
of some observable algebra, then the scaling behavior of A(L b ) may show up in a 
positive scalar factor Q, 

A(s ■ x) = s K ■ A(x) 
A(x) = Q(x) ■ A'{x) 

Q(x) cx x K (86) 

where the algebra element A' docs not show significant scaling behavior, or it's 
average might even vanish due to symmetry. The following generalized multifractal 
analysis applies to such observables. With this restriction in mind we proceed as 
follows. 



4.1. Normalization Exponents 

To study the statistics of box-observables one takes a large number N of equivalcntly 
prepared systems and records the values of the box observable. This leads to the 
set {Qi{Lb)} i=l ft- The mean value of a function of Q is then calculated as 

1 * 

i 

A corresponding smooth distribution function n(Q,L^, N) is assumed to describe 
average values (provided TV is large enough), 

oo 

(F(Q(L b ))) ft = JdQ F(Q)U(Q, L b , N) . (88) 
o 

In the absence of intrinsic length scales, i.e. compared with L b intrinsic length 
scales are either much less or much larger than Lb, one expects power law scaling 
to occur in averages. To make contact to the multifractal analysis of measures at 
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an early state of the investigations let us introduce a symbolic system size L by 
reparametrizing the number N as 



L:=[LZN) (89) 

Let N < N be the number of non- vanishing values for Qi(L b ), define the normal- 
ization of Q by 

N 

\\Q(L b ,L)\\:=YlQi( L b) (90) 

i 

and, by keeping L fixed in the averaging procedure, identify the normalization 
exponent XW in the regime of power law scaling via 

||Q(L 6 ,L)||ocLf WI . (91) 

X^ can be any real number. The crucial step is the construction of a normalized 
box- observable P([Q], L) associated with Q: 



Pi([Q],L b ,L)= n^^H . E«i*£)=l. (92) 

Pj([Q], Lfc, L) is normalized with respect to the symbolic system of linear size L but, 
in general, not additive, 

N 

P{[Q],L)^^Pi{[Q},L h ,L). (93) 

i 

In the usual multifractal analysis one imagines the box-observables to be measures 
for which the normalization condition is very natural. Nevertheless, when proving 
the general features of the functions r(q) and f(a) the normalization condition was 
essential, but no use was made of additivity. By fixing L we are able to attach to 
the positive, normalized box-observable P([Q], L b , L) the multifractal analysis by 
defining r^(q), a^(q) and f^(a^) with the help of the partition sum, 

N 

Z([Q];q,\):=J2P?([QlL b ,L), (94) 

i 

in direct analogy to Eqs. (|58|). For these functions all the general features explained 
in Sec. 3 follow if the ratio A = L b /L is choosen small enough, and the observables 
Qi(L b ) are in the scaling regime. The results of Sec. 3 can be translated to the 
following statements 

• For positive box-observables Q{L b ) the scaling behavior of its distribution 
function can be described by the corresponding normalization exponent X^ 
and a single-humped function f^(a^') 



N 
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with 

£3= (96) 
In A 

The scaling relations for the moments and the typical value, 

Qtyp(L b ) := exp{(lnQ(L b ))} , (97) 

read 

(P«([Q],L b ,L)) oc Lf+- tQ1 W 

(Q q (L b )) oc L* W , x(q) =D + TM(q) + q-XM 

Q typ (L 6 ) oc L™° +X (98) 

The statements of Sec. 3 about properties of r, a, / such as monotonicity, 
curvature, Legendre transform, etc. do also apply for T^\a^ and In 
general one expects D — d, unless a systematic scaling dependence of N(N) 
appears. 



4.2. Local observables 

In case that Q is a local additive observable, i.e. for each subdivision of a starting 
system of linear size L into N boxes of linear size L b , Q(L b ) is an observable and 
additivity holds, 

N 

Q(L) =J2Qi(Lb), (99) 

i 

the generalized multifractal analysis can also be applied to the investigation of 
scaling properties with respect to the system size L which is then a true system 
size. We focus on the situation where L b and L are in a regime of power law scaling 
and vary both lengths independently. We thus have to be aware of an additional 
normalization exponent with respect to L which is, in general, different from 

\\Q(L b ,L)\\o,Lf Q] L- YlQ] . (100) 

The box-probabilities studied in Sec. 3 fall, of course, into the class of local observ- 
ables with trivial normalization exponents X = Y = 0. For local observables one 



finds, besides the scaling relations of Eqs. (95, Bq) with respect to L b , (by varying 



L and taking the normalization exponent Y^ of Eq. ( |100| ) into account) 

mLb )) L oc , Q typ (L b ,L) oc z^'VK^wi) (101) 

with 

x(q) =D + rW(q)+qXW , y(q) = D + (q) + qY^ . (102) 

As an illustration let's discuss two examples. A trivial example is the box-observable 
Q := PLff L~ Y where P is a box-probability with multifractal functions t, a and 
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/ and X,Y are real numbers. Then, = r, cJ Q J = a, = /, and = X, 
y[Q] = y < As another example consider the box-observable Q = P m where P is 
a box-probability and m some real number. All the relevant scaling information is 
already contained in the r(g) function of P. Thus, it is easy to check that 

X [P™\ =Y [P™\ = r(m) 

o^ pm \q) = m ■ a(q • m) - 
/[ pm ](a [pm] ( 9 )) = (gm) • a(gm) - r( ? m) . (103) 

The last of these three equations yields, due to the Legendre transform property of 
f{a) with respect to r{q), 

(aM) = f(a(qm)). (104) 

It is tempting to expect from this examples that /(a)-spectra of different observ- 
ables, defined for one ensemble of equivalently prepared systems, are either identical, 
if only the normalization is distinct, or are related in a simple manner, if powers of 
elementary observables are involved. 

4.3. Typical Scaling Variables 

In the vicinity of a phase transition one expects a universal behavior of the scal- 
ing exponents. They may depend on general properties of the systems dynamic, 
however not on microscopic details. We have already observed that different wave 
functions in the LD transition regime showed the same /(a)-spectrum. A further 
indication for universality of /(a)-spectra was provided by the observation of Huck- 
estein and Schweitzer 43 that the local equilibrium current density as well as the local 
magnetization 52 show the same /(a)-spectrum as the wave functions in a quantum 
Hall system. However, universality of scaling variables in the LD transition (which 
are non equilibrium properties) is even more interesting. 

Recall the phenomenological scaling theory of Abrahams et al. for the LD tran- 
sition (see Sec. 2.4) where the conductance g is assumed to be a scaling variable 
on all length scales. As we already mentioned in Sec. 2.4 the conductance of meso- 
scopic systems has a broad distribution which cannot be described in terms of the 
mean value of g. Consequently, one has to consider the distribution function TL(g, L) 
which can be characterized by relevant parameters g rc i for which (3- functions exist. 
If g is choosen for the f(a) analysis a candidate for g re \ is obviously given by the 
typical value gt yP (L) = exp {(lng(L))} 

gty P (L)-g; yp KL al ° 9] + xlsl (105) 
and by comparing Eqs. ( ^7|j98| ) one concludes 

l/u = c$+XW. (106) 



The universality of v suggests the /(a)-spectra of normalized scaling observables be- 
ing universal, too. Notice that Eq. (106) provides a method to calculate the critical 
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Fig. 4. Thouless numbers for 6 realizations of a quantum Hall system with linear size 120r c 
restricted to the lowest Landau level. LU C is the cyclotron frequency. Energy windows for data 
evaluation as described in the text are indicated by vertical lines (after Ref. 23). 

exponent v with only the help of multifractal critical numbers. This gives the possi- 
bility for calculating v by a different method than the finite size scaling method. 14,15 
In addition, if both methods coincide the multifractal analysis demonstrates that 
typical observables are relevant for scaling; a consequence of the fluctuating self 
similar structure underlying multifractal objects. 

In order to check the predictions of the multifractal analysis Fastenrath et al. 23 
calculated numerically Thouless numbers for the QHS and applied the multifractal 
analysis to the data. The Thouless numbers are defined by the ratio of the energy 
shift due to a change of boundary conditions (e.g. periodic — ► anti periodic) with 
respect to the mean level spacing. They can be thought of as being a transport 
quantity showing the same qualitative (and scaling) behavior of a (dimensionless) 
dc conductance. 53 

In Fig. ^ data for the Thouless numbers denoted by g are shown which are 
calculated for the lowest Landau level of 6 realizations of a QHS with linear size 
120r c . One can see the critical regime outside of which the states are localized and 
have zero dc conductance. In the critical regime there are large fluctuations and the 
corresponding histogram of Fig. |^ demonstrates the broadness of the distribution 
function. To apply the multifractal analysis, Thouless numbers g{L) for a large 
number of different systems with varying system sizes from 30r c to 120r c were calcu- 
lated. As the box observable Q the modulus of the deviation Q(L) := \g(L) — g* yp \ 
was chosen where g^ yp was indeed size independent. Fig. [| shows that the cor- 
responding /(a)-spectrum of the Thouless numbers is very similar to that of the 
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Fig. 5. Histogram of Thouless numbers calculated for the lowest Landau level of one quantum 
Hall system. 

box-probabilities. 20 Data have been taken from a narrow window in the critical 
region as indicated in Fig. ^. The authors find 

aff 1 = 2.22 ± 0.05 , X [9] = -1.75 ± 0.05 , v = 2.2 ± 0.3 . (107) 

They also calculated v directly from | lngtyp — ln<7 t yp*| oc L 1 ^ and found 

z/ = 2.3±0.2. (108) 

In either case they used g^ yp — 0.27 ± 0.03, determined at the band center. The 
data for v are in agreement with the high precision value of Huckestein and Kramer, 
obtained by finite size scaling calculations. 51 

The multifractal analysis of Fastenrath et al. 23 confirm the one-parameter scaling 
theory of the LD transition in terms of the typical conductance and demonstrate 
the broadness of the conductance distribution at criticality. 

4.4. Bounds for the Correlation Length Exponent 

In this section we show how the multifractal analysis provides a lower bound for the 
critical exponent v. We refer to the LD transition for concreteness but the results 
generality is reflected by the fact that only a few assumptions such as one-parameter 
scaling are needed. To obtain upper bounds for v more restrictive assumptions are 
needed. 

That v has a lower bound given by 

2 

v > - (109) 
a 
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Fig. 6. /(o) spectra of normalized (right) and unnormalized (left) Thouless numbers. The 
displacement of the cxtrcma is just the normalization exponent X . The broken line to the right 
is a parabolic approximation (after Rcf. 23). 



where d is the Euclidean space dimensionality is known since the work of Chayes et 
al. , 54 A trivial upper bound for v is given by 

v < oo (110) 

otherwise the definition of the exponent would be meaningless (e.g. in the Kosterlitz- 
Thouless transition). 55 We give arguments which improve both bounds. Our argu- 
ments rely on the assumption that one parameter scaling holds true in the vicinity 
of the LD transition point, on the analytic properties of r{q) functions and on the 
universality hypothesis. 

In a field theoretical statistical model for a critical phenomenon as outlined in 
Sec. 2 a scaling operator couples to the scaling field t (in the LD transition 
problem t corresponds to \E — E c \ =: t) in the Hamiltonian 




(111) 



Here S(r) is a local scaling operator of the field theoretical statistical model. The 
scaling dimensions y(n) corresponding to the scaling operators 




(112) 



are defined with respect to the renormalization group action (see e.g. Ref. 56). 
Renormalizability of the theory then means that only a finite number of these 
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exponents are positive (relevant) and most of them are negative (irrelevant). Here 
we identify them via the finite size scaling properties of the average 

^Jd d rS n (r)^J oci K(n) (113) 

where the statistical model is defined on a finite system fi of linear size L. For 
a relevant scaling operator the integral diverges for L — > oo and consequently the 
scaling fields tW in tM • have to be zero, i.e. the critical point is reached. 
One-parameter scaling means that there is exactly one relevant operator and thus 

«(n) < for n > 2. (114) 



If by accident or symmetry properties the expectation value in Eq. (113) vanishes 
we assume that there exists a corresponding positive scalar observable Q which has 
the same scaling behavior as S(r), in the spirit of Eq. (|§6|) . If this assumption holds 
true, then the general analytic properties of r(q) functions tell that 



dn(n) 
dn 
d 2 n(n) 



dr* ^ ■ < 115 > 

Consequently the scaling dimensions have positive curvature. Since the exponent 
v is given by y(l) = \jv (for consistency with the divergence of the correlation 
length) we have 

X [B] = -l/v (116) 

In order to have n{n) monotonically decreasing guaranteeing renormalizability the 
lower bound 

v > 1/DW (in) 

has to be respected. Furthermore, the one-parameter scaling condition Eq. (|ll4j) 
requires 

V> 1MW- (118) 



This resembles and improves the Chayes et al. criterion Eq. (10£). The observation 
of universality in f(a) spectra for different observables suggests D^\q) = D(q). 
Relying on the universality hypothesis for /(a) the multifractal analysis of the wave 
function in the LD transition problem already allows to obtain a lower bound for 
v, v > 2/D(2). 

In Sec. 4.3 we have found that it is the typical conductance gt yp rather than the 
average conductance (g) which serves as a one-parameter scaling variable for the 
LD transition and scales like 

\g typ (L) - g; yp \ oc L l ' v . (119) 
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The corresponding /3-function has an unstable fixed point. The exponents x i9 ^{q) 
describing the scaling of moments (\g(L) — fftypl 9 ) K L xl9](q ^ are given by 

x^(q) = D + T^(q)+qX^ . (120) 

The analytic properties of x^ are summarized as follows. 

= aW(q)+XW 

dq 

d 2 x^(q) da^(q) 



dq 2 dq 



<0. (121) 



Since a;[ ff l(0) = and 1/v = affl + X^ > the exponents x^(q) are positive for 
small positive values of q. For large enough values of q they will become negative if 
D\£+XW <0 which means that high moments of \g(L) — g% yp \ may decrease when 
approaching the transition point. This is not in contradiction to renormalizability 
and does not yield a rigorous upper bound for v. However, if at least the first 
moment has positive scaling exponent, x^(l) = D + X^ > 0, then an upper 
bound for v can be concluded 

v<-r^- . (122) 

af-D 

Note, that this upper bound relies on the additional, though intuitive, assumption 
about a positive scaling exponent for the mean deviation (g(L) — .9t yp )- 

We give three examples for which the validity of the bounds can be checked 
relying on the universality of f(a) spectra. 

1. For the QHS the values for v, D(q) and «o reported above are compatible with 
the bounds, 

1.3 < v = 2.3 < 3.3. (123) 

2. A similar conclusion holds for the 2-d spin-orbit model (the calculation of v = 
2.75 was done by finite size scaling methods 11 and a = 2.175 was calculated by 
the multifractal analysis 57 ), 

1.2 < v = 2.75 < 6 (124) 

3. The one-loop results of Wcgner's non linear er-model 47 in d = 2 + e leads to a 
parabolic /(a)-spectrum as shown in Ref. 19. The results are 

v = lie, a = d + e, D{2) = d- 2e + (125) 

Thus, the one loop result meets the upper bound and the lower bound is valid up 
to e < 2/3. Higher loop orders 47 lead to f(a) spectra which violate the condition 
of constant curvature if taken seriously for arbitrary values of e < 1. This indicates 
that the higher orders are improvements only for very small values of e (as often 
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happens in asymptotic expansions). The applicability of the loop expansion is an 
open question. 

5. Correlations in Multifractals 

In this section we ask for correlations in positive, local observables as introduced in 
Sec. 4.2. For example, a local field <p(r) gives rise to local box-observables 

Qi{L b ) = J d d r\ <p(f) | . (126) 

box(i) 

We are interested in the scaling properties of 

M^(r,L b ,L) := <Qf (L b )Q q +s (L b )) L (127) 
where the average is over all pairs of boxes with fixed distance r = sL b . 

5.1. Scaling Relations 

Usually, in critical phenomena one studies correlations for infinite system size (and 
L b being microscopic) as a function of r alone. Here we consider a regime, where 
both r and L are able to show scaling behavior. We thus start with an ansatz 

M [q \r,L b ,L)ozL x b 2iq) L- y2iq) r~ z{q) (128) 

and try to find relations of the exponents X2{q), 2/2(9), z(q) to the previous ones 
x(q) and y(q)l A typical heuristic scaling argument to fix scaling relations relies 
on considering limiting cases where scaling will be already violated. However, since 
the limit is reached continuously and scaling exponents are unique one can conclude 
scaling relations due to consistency. We thus proceed by considering the limiting 
cases for M™: (i) r — L b and (ii) r — L. In case (i) we expect the asymptotics 

oc (Q 2q (L b )) L oc Lf^L-vi^) . (129) 

In case (ii) we expect the asymptotics 

cx (Q^(L b )Q^(L b )) L cx {Q q (L b )) 2 L oc L^L"^') (130) 

where Q q (L b ) is the value for Q q (L b ) already averaged over one system of size L. 



Comparison with Eq. (128) yields the scaling relations 



V2{q) = y(2q) 
x 2 {q) = x(2q) 
z(q) = 2y(q) - y(2q) = 2x(q) - x(2q) (131) 

which coincide with those of a more sophisticated derivation by Cates and Deutsch 58 
in the case of Q being a local measure. Notice that the normalization exponents X 

tin the following we suppress the superscript [Q] in exponents. 
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and Y cancel in z(q) leading to z(q) = D + 2r(q) — r(2q). The analytic properties 
of z{q) 1 X2{q) and 1/2(1) are summarized as follows. z(q) vanishes at q = where 
it has a minimum. It is monotonically increasing (decreasing) for q > (q < 0), 
and it's slope vanishes for q — > ±00. £2(9) and 1/2(1) have negative curvature. 
Their monotonicity depend on the sign of a(2q) + X and of a(2q) + Y, respectively. 
Since a (q) is monotonically decreasing they are monotonically increasing provided 
Doe > —X, — Y. Notice that X2 (<?) and 2/2(9) depend on the normalization exponents 
and may become negative for positive values of q. 

The fact that can depend on the system size L is unusual compared to 

ordinary critical phenomena where the correlation function, the susceptibility x( r )i 
is expected to depend on the distance r only, provided the thermodynamic limit is 
reached and the correlation length £ is infinite at the critical point T c . We like to 
comment on this in the following. 

In Sec. 3.2 we discussed an interpretation of multifractality of the box-probability 
in the context of the LD transition: Multifractality reflects broadness of the box- 
probabilities distribution function on all length scales which is due to the dependence 
of the box-probability in each box on a large number of conditions, simultaneously. 
More general, we call a situation, where a local box observable depends on a large 
number of conditions for the entire system of linear size L, simultaneously, a situ- 
ation of "many parameter (MP) coherence" . In the context of the LD transition 
coherence at zero temperatures is due to quantum mechanical phase coherence of 
the electrons wave function and disorder introduces a huge number of parameters, 
e.g. the position of point-scatterers. It may happen that MP coherence is valid 
only up to a certain scale L < L, defined implicitly by being independent of L 
for L > L. We call L a MP coherence length. If such crossover in M™ exists, two 
situations have to be distinguished. First, L introduces a cut-off for correlations. 
For example, the correlation length £ is a MP coherence length of this kind. Al- 
ternatively, L does not introduce a cut-off and correlations still show homogeneity 
exponents for distances r > L. An example for this kind of MP coherence length 
in the LD transition problem will be discussed below. In the latter situation the 
following scaling behavior of is expected to occur for r <C L <C L, 

AfW (r, L, L) oc L-f=(9) r -*W , (132) 

and for L ^ r 3> L, 

M^{r,L,L) ocr^, (133) 

respectively. Thus, for r 3> L the usual behavior is recovered. However, for this sit- 
uation the multifractality on scales less than L is still reflected by the (/-dependence 
of z(q) the scaling relation of which can be concluded by a similar reasoning as that 



leading to Eqs. (131) 



S(q) = y 2 (q) + z(q) = 2y(q) . (134) 

The analytic properties of z(q) are such: z(q) has negative curvature and its mono- 
tonicity properties depend on the sign of ot(q) + Y. It can thus happen, that z(q) 
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is negative for a wide range of q values. A system which is never MP coherent, 
i.e. L is microscopic, the single-fractal situation applies and, then, z(q) = and 
~z(q) = 2y(q)=y(2q). 

5.2. Application to LD Transitions 

In the LD transition problem the density of states of an individual system is a 
suitable candidate for a local box-observable. It is defined by 

p(E, r) := (r\S(E - H)\r) = ]T \Mr)\ 2 S(E - e a ) (135) 

a 

where tpa are eigenstates of H with respect to energy e a . One peculiar feature of 
the LD transition is that the average density of states, p{E) = -^a-Tr 8{E — H), is 
not an order parameter and is L-independcnt. The corresponding box observable 
is 

Q P (L b ,E):= J d d rp(E,f) (136) 

box(i) 

The scale-independence of p{E) determines the normalization exponents to be 

X p = 0, Y p = -D. (137) 
The fact that the average of p(E) is not an order parameter (/3 = 0) is equivalent 

y 

'typ 



to x p (l) = 0. However, the typical value of Qt yp can serve as an order parameter, 



because 

Q p typ cx L-(«°-^ (138) 
and consequently, by L approaching £ c , we have 

Q p typ oc t u < a °- D ^ . (139) 

Thus, in contrast to the mean value (Q p ), Qt yp 1S a ble to indicate the LD transition 
and we have reason to call 

Ay P ~ K^o + Y) = v{a - D) (140) 

the typical order parameter exponent. 

The exponents z(q) do not depend on the normalization exponents and are given 
as z p (q) = D + 2r(q) — r(2q) whereas the exponents z p (q) do depend on Y p = —D, 

~z p {q) = 2(q-l){D{q)-D) (141) 

have negative curvature, are positive for < q < 1, vanish at q = 0, 1 and are 
negative elsewhere. The function r(q) is, due to Eq. ( |135| ), the same as for the wave 
function itself. 

For q — 1 we can compare with the result of Chalker and Danicll 30 for the scaling 
of the density correlator at the LD transition of a quantum Hall system, 

{p(E + uj/2,f)p(E-uj/2,r')) , (142) 
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where u) sets the scale of a MP coherence length = (mp)^ 1 ^ 2 . L u can be 
interpreted as the linear system size for a system with mean level spacing of about 
uj. 30 Their result is: for r <C L w correlations scale like r~ 38 whereas for r 3> L u 

The first of these results was already discussed in 



Sees. 2,4 and can be recovered from Eq. (131) 



z = z p {\) = D -D(2) = 0.38. (143) 

The second result is consistent with the interpretation of z being the distance ex- 
ponent in a regime r 3> L u , i.e. 

z = z p (l)=0. (144) 

A similar observation was made in Ref. 12 for the LD transition in 2-d system with 
spin-orbit coupling , though the results were not conclusive concerning quantitative 
results for z p (l). Furthermore, the scaling of the density correlator with respect to 



L^, fulfills the scaling relation of Eq. (131) which means here 



l£(l)=y"(2) = -z = -0.38. (145) 

This came out in the work of Chalker et al. by studying the combined variable L w k 
where k is the inverse wavelength in the Fourier-transform of the density correlator 
and probes long distances in the limit k — » 0. Instead of interpreting ui via L u one 



can (more directly) look at the density correlator of Eq. (142) as describing the 
density correlations in energy and, owing to Fourier transformation for the limit 
to — > 0, the long time correlations of wave packets. Huckestein and Schweitzer 50 
checked the latter point of view by explicitly calculating the long time correlations 
of wave packets in the critical regime of a quantum Hall system and found excellent 
agreement, i.e. z = 0.38 ± 0.02. 

Unfortunately, no results for values q ^ 1 for z p (q), y^iq) and z p (q) are available 
up to now. Thus, the picture developed here needs further tests. 

5.3. Conformal Mapping in 2-d 

We apply conformal mapping arguments to find a relation between finite size scaling 
(FSS) methods relying on strip- like systems and the multifractal analysis for square- 
like systems in 2-d. To begin with some ideas about conformal mapping are reviewed 
(for a satisfactory treatment see e.g. Ref. 59). 

For critical correlation functions A(r) oc r~ K the following assumption seems 
plausible. Scale invariance, as reflected by Eq. ([!]), should hold also for local trans- 
formations which preserves angles but may change scales locally in the sense that 
correlations of one system with a certain geometry are mapped onto those correla- 
tions of a similar system the geometry of which is determined by the transformation. 
Such transformations are called conformal mappings. The corresponding Jacobian 
J has to fulfill 

v ■ w Jv ■ Jw 

m = T75 M 6 

(iPvP) 1/2 ((Jv) 2 {Jw) 2 ) 1/2 
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where v, w are vectors of the tangent space at a given point. In 2-d any holomorphic 
function f(z = x + iy) is a conformal mapping. Especially, the complex logarithm 
provides a conformal mapping of the entire plane onto an infinite strip with periodic 
boundary conditions. Let's introduce Cartesian (x,y) and polar coordinates (r, tp) 
on the plane by x + iy = re ltp . Then the conformal mapping F(x + iy) — ^- \n{x + 
iy) =: u + iv maps a ring area, denoted by C, with inner radius l/L and outer 
radius L onto a strip of width M with length V = — InL, centered at u = and 
fulfilling periodic boundary conditions in ^-direction. For a correlation function on 
the plane (x,y), A{r) oc r~ K which behaves regular for r — > we can express the 
exponent k by an integral expression in coordinate free terms by (cf. Ref. 60) 



-1 
2^ 



lo (147) 



c 



where the two- form uj := d * d(ln A), * is the Hodge-operator corresponding to the 
Euclidean metric, and l/L is assumed to be small enough to give vanishing contri- 
bution to k by applying Stokes' theorem. By lifting uj with the conformal mapping 
F, and assuming exponential decay of the corresponding correlation function along 
the strip, A(u, v) oc exp —\u\/£(M), one arrives at the result 

^Ml = ±. (148) 

We wish to apply this result to correlation functions of multifractal correlations 
AfW with L b being microscopic. However, in the regime r <C L, L conformal map- 
ping arguments of the kind presented here can never apply since a second length 
scale besides the distance r appears. We thus have to focus on the regime r ^S> L 
where behaves as 

AfW(r)ocr -s ^. (149) 

At the critical point of the LD transition, where power law in the plane is valid, the 
decay lengths ^ (M) of the corresponding correlations in the strip-geometry vary 
linear with M. This fact is exploited in the FSS analysis of the LD transition (see 
e.g. Ref. 14). Thus, the FSS variable A^(M) := 2fM (M)/M*becomes a constant 
A[?' at the critical point. 

Making the hypothesis, that the conformal mapping result Eq. (148) applies 
to the multifractal correlations in the regime described by Eq. ( |149 ) one would 
conclude 

AM = 2/(tt%)) . (150) 

However, in the situation of the LD transition z{q) takes negative values for q > 1 
which is counterintuitive and suggests exponential growth of correlations in the strip 
geometry. Furthermore, the exponents z(q) and the decay lengths ^ are attached 
to averages of powers of the correlation function, though the precise statement on 



tThe factor of 2 is due to convention. 
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the existence of the decay length in the strip-geometry relies on the average of 
the logarithms of the correlations (see e.g. Ref. 14). Thus, we leave the validity 
of Eq. (15C) as an open question and turn over to the typical values defined by 
the average of logarithms. In FSS calculations it is the average of the logarithm 
of the correlation function which is calculated by knowing that this quantity is 
self-averaging. Therefore, we denote results of such FSS calculations by A* yp and 
a conjecture about the relation to multifractality (relying on conformal mapping 
arguments) reads 

A^p = 2/(7rz typ ) , ~z typ = 2{a + Y) (151) 

where the expression for z typ is a consequence of Eqs. (134,101). 

In the context of the LD transition we found Y p = — D for the local observable 
being the density of states. This observable involves the square amplitude of wave 
functions and localization lengths are usually defined with respect to the modulus 
of the wave function. This explains the conventional factor of 2 in the definition of 
A(M). 

Up to now, only data for two universality classes are known where one can check 
the prediction of Eq. (151): (i) quantum Hall systems and (ii) 2-d Anderson model 
with spin-orbit scattering. Taking into account that present days calculations, lack- 
ing larger systems, are able to produce results for Aj: yp and ao — D with a precision 
hardly exceeding 5%, the results are not in conflict with Eq. ( |151[ ): (i) A* yp = 1.14 61 
and l/(7r(oio - 2)) = 1.06. 43 (ii) A typ = 1.88 62 and l/(7r(a - 2)) = 1.82. 57 



6. Conclusions 

After reminding some aspects of critical phenomena in Sec. 2 we described the 
multifractal analysis of broadly distributed observables in the critical regime of a 
critical phenomenon. Broadness of the distribution on length scales much less than 
the correlation length £ c is determined by the single- humped /(a) spectrum of 
normalized box-observables. The parabolic approximation for /(a) corresponds to 
a log-normal distribution. Scaling of the q-th moment is given by the function r(q) 
defining generalized dimensions D(q) = T(q)/(q — 1) (D = D(0) being the geometric 
fractal dimension) and by the normalization exponents AT, Y . The function r(g) 
is monotonically increasing and has negative curvature. It is related to /(a) by 
Legendre transformation (Sees. 3 — 4). Since the distribution is broad it is not 
possible to characterize it by the mean value. Instead, the typical value defined as 
the geometric mean is a self-averaging quantity which scales with exponent ao + 
X where ao is the maximum position of fipt), /(ao) = D. We interpreted the 
appearence of broad distributions being due to many parameter coherence, saying 
that local observables depend on a large number of conditions for the entire system, 
simultaneously. 

Furthermore, correlation functions of broadly distributed observables at critical- 
lity show scaling dependence with respect to the distance between local observables 
and with respect to the linear system size L or, equivalently, with respect to a length 
L indicating the range of many parameter coherence. For distances r exceeding the 
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length L correlations will only depend on r (Sec. 5). Although compatible with 
known results, this scenario needs further investigations. 

Scaling exponents of correlations are related to r(q), X, Y by scaling rela- 
tions (Eqs. (131 134)) Therefore, a complete characterization of critical exponents 
in terms of multifractal spectra is possible. 

We applied the multifractal analysis to the localization-delocalization transition 
induced by disorder (LD transition) leading to the following conclusions. The formal 
order parameter, the average density of states, does not show the LD transition 
which means the order parameter exponent vanishes, (3 — 0. Consequently, usual 
scaling relations tell that correlations show a distance exponent z = 0. In contrast 
to the average density of states, the typical value is able to reflect the LD transition 
and the corresponding typical order parameter exponent and distance exponent 
are /?ty P = (c*o ~ D)v and Zt yp = 2(ao — D), respectively (Sec. 5). Here v is the 
critical exponent of the localization length. By choosing the typical conductance 
as a scaling variable an expression for v in terms of only multifractal exponents 
can be given, v = («o + X^ 1 ')^ 1 . The multifractal analysis yields lower as well as 



upper bounds for v in terms of D(2) and ao, respectively (Eqs. (118, |122| )). The 
exponent z describing the distance exponent in the regime r <C L, L is given by 
z = d — D{2) ^ z which is excellently confirmed for the LD transition in quantum 
Hall systems (Sec. 5). 

With the help of conformal mapping arguments in d — 2 we suggested a relation 
between the critical value of the (typical) renormalized localization length A c of 
strip-like systems (being the scaling variable in finite size scaling calculations) and 
the exponent a$ (Eq. (|l5l| )). This relation is in agreement with presently avail- 
able data. A corresponding relation between moments (Eq. ( |l50|) ) needs further 
numerical investigations. Such relation is highly desirable since it would allow to 
demonstrate complete equivalence between the multifractal analysis and the finite 
size scaling approach. 

Finally, we comment on the question of how to determine the critical point T c 
by the multifractal analysis. For finite system sizes L the critical regime of values T 
around T c , characterized by L <C £ c , is finite. Within this regime the /(a) spectra 
can be calculated by determining linear regimes in log- log plots (see Eqs. (|58|)), 
These values will slightly differ from the universal values at T c . The difference 
(finite size corrections) will scale, e.g. 63 a (T) - a (T c ) oc \T - T C \L X I V . Thus, 
in principle, one can determine T c and \Jv from finite size corrections. 17 There 
have been speculations 64,65 to determine T c , in the context of the LD transition, by 
requiring a certain fractal dimension to coincide with the lower critical dimension 
di which is believed to be di — 2. Since it is known that there is a spectrum of 
fractal dimensions, such criteria are not evident. To the authors knowledge the only 
motivation for conjecturing -D(l) = di as a criterium to fix T c comes from the one- 
loop result of Wegners non-linear er-model in d = 2 + e 19 ' 47 where D(q) = 2 + e — qe. 
However, this approximation is expected to deviate from exact values as e — > 1. 
Thus, in order to determine T c by the multifractal analysis one has to establish 
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universality in multifractal exponents rather than to establish certain values for 
these exponents. 

Acknowledgements 

This work was supported by the Deutsche Forschungsgcmcinschaft. I thank the 
Department of Theoretical Physics in Oxford, U.K., for the kind hospitality during 
a visit in which most part of this work was performed. I gratefully acknowledge 
helpful discussions with J. Chalker, U. Fastenrath, J. Hajdu, B. Huckestein, W. 
Pook, L. Schweitzer and M. Zirnbauer. 

References 



1. P. W. Anderson, Phys. Rev. 109, (1958) 1492. 

2. P. A. Lee, T. V. Ramakrishnan Rev. Mod. Phys. 57, (1985) 287. 

3. A. J. McKane, M. Stone, Ann. Phys. 131, (1981) 36. 

4. R. Johnston, H. Kunz, J. Phys. C16, (1983) 3913. 

5. E. Abrahams, P. W. Anderson, D. C. Licciardello, T. V. Ramakrishnan, Phys. Rev. 
Lett. 42, (1979). 

6. J. T. Chalker, M. Bernhardt, Phys. Rev. Lett. 70, (1993) 982. 

7. U. Fastenrath, Solid State Commun. 76, (1990) 855. 

8. K. von Klitzing, G. Dorda, M. Pepper, Phys. Rev. Lett. 45, (1980) 494. 

9. R. Prange, S. Girvin (eds.), The Quantum Hall Effect, (Springer, New York, 1990). 

10. S. N. Evangelou, T. Ziman, J. Phys. C20, (1987) L235. 

11. U. Fastenrath, G. Adams, R. Bundschuh, T. Hermes, B. Raab, I. Schlosser, T. Wehner, 
T. Wichman, Phystca A172, (1991) 302. 

12. J. T. Chalker, G. J. Daniell, S. N. Evangelou, I. H. Nahm, J. Phys.: Condens. Matter 
5, (1993) 485. 

13. U. Fastenrath, Phystca A189, (1992) 27. 

14. J. L. Pichard, G. Sarma, J. Phys. C14, (1981) L127. 

15. A. MacKinnon, B. Kramer Phys. Rev. Lett. 47, (1981) 1546. 

16. B. Kramer, K. Broderix, A. MacKinnon, M. Schreiber, Physica A167, (1990) 163. 

17. J. Bauer, T. M. Chang, J. L. Skinner, Phys. Rev. B42, (1990) 8121. 

18. M. Schreiber, H. Grufibach, Phys. Rev. Lett. 67 (1991) 607. 

19. C. Castellani, L. Peliti, J. Phys. A19, (1986) L429. 

20. W. Pook, M. Janfien, Z. Phys. B82, (1991) 295; W. Pook, PhD thesis, Universitat 
zu Koln, 1991. 

21. B. Huckestein, L. Schweitzer, in High Magnetic Ftelds in Semiconductor Physics 
HI: Proceedings of the International Conference, Wiirzburg 1990, ed. G. Landwehr 
(Springer Series in Solid State Sciences 101, 1992) p. 84. 

22. B. Huckestein, B. Kramer, L. Schweitzer, Surface Science 263, (1992) 125. 

23. U. Fastenrath, M. JanBen, W. Pook, Physica A191, (1992) 40. 

24. A. W. W. Ludwig, Nucl. Phys. B330, (1990) 639. 

25. B. Duplantier, A. W. W. Ludwig, Phys. Rev. Lett. 66, (1991) 247. 

26. K. B. Efetov, Adv. Phys. 32, (1983) 53. 

27. A. M. M. Pruisken, in Ref. 9, p. 119. 

28. F. Wegner, Z. Phys. B51, (1983) 279. 

29. H. E. Stanley, Phase Transitions and Critical Phenomena (Clarendon, Oxford, 1971). 

30. J. T. Chalker, G. J. Daniell, Phys. Rev. Lett. 61, (1988) 593. 



Multifractal Analysis ... 39 



31. J. T. Chalker, J. Phys. C20, (1987) L493. 

32. B. L. Al'tshuler, V. E. Kravtsov, I. V. Lerner, Sov. Phys. JETP 64, (1986) 1352. 

33. B. L. Al'tshuler, V. E. Kravtsov, I. V. Lerner Phys. Lett. A134, (1989) 488. 

34. B. Shapiro, Philos. Mag. B56, (1987) 1031. 

35. B. B. Mandelbrot, The Fractal Geometry of Nature (W. H. Freeman, San Francisco, 
1983). 

36. U. Frisch, G. Parisi, in Turbulence and Predictability of Geophysical Fluid Dynamics 
and Climate Dynamics, Proceedings of the International School of Physics "Enrico 
Fermi" LXXXVIII, Varena, ed. M. Ghil (North Holland, New York, 1985). 

37. P. Meakin, A. Coniglio, H. E. Stanley, T. A. Witten, Phys. Rev. A34, (1986) 3325; 
M. E. Cates, T. A. Witten, Phys. Rev. A35, (1987) 1809. 

38. T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, B. I. Shraiman Phys. Rev. 
A33 (1986) 1141. 

39. B. B. Mandelbrot, J. Fluid Mech. 62, (1974) 331. 

40. H. G. E. Hentschel, I. Procaccia, Physica D8, (1983) 435. 

41. T. Tel, Z. Naturforsch. 43a, (1988) 1154. 

42. A. Chhabra, R. V. Jensen, Phys. Rev. Lett. 62, (1989) 1327. 

43. B. Huckestein, L. Schweitzer, Helv. Phys. Acta 65, (1992) 317. 

44. J. T. Chalker, J. Phys. C21, (1988) L119. 

45. H. Aoki, J. Phys. C16, (1983) L205. 

46. H. Aoki, Phys. Rev. B33, (1986) 7310. 

47. F. Wegner, Z. Phys. B36, (1980) 209; F. Wegner, Nucl. Phys. B316, (1989) 663. 

48. G. Paladin, A. Vulpiani, Phys. Rep. 156, (1987) 147. 

49. S. Hikami, Prog. Theor. Phys. 76, (1986) 1210. 

50. B. Huckestein, L. Schweitzer, Phys. Rev. Lett. 72 (1994) 713. 

51. B. Huckestein, B. Kramer, Phys. Rev. Lett. 64, (1990) 1437. 

52. L. Schweitzer, 1992 (private communication). 

53. D. J. Thouless, Phys. Rep. 13C, (1974) 93. 

54. J. T. Chayes, L. Chayes, D. S. Fisher, T. Spencer, Phys. Rev. Lett. 57, (1986) 299. 

55. J. M. Kosterlitz, D. J. Thouless, J. Phys. C6, (1973) 1181. 

56. F. Wegner, in Phase Transitions and Critical Phenomena, Vol. 6, eds. C. Domb, M. 
S. Green, (Academic, London, 1976). 

57. L. Schweitzer, 1993 (private communication). 

58. M. E. Cates, J. M. Deutsch, Phys. Rev. A35, (1987) 4907. 

59. J. L. Cardy, J. Phys. A17, (1984) L385. 

60. J. L. Pichard, G. Sarma, J. Phys. C14, (1981) L617. 

61. B. Huckestein, 1993 (preprint). 

62. U. Fastenrath, 1992 (unpublished). 

63. B. Huckestein, L. Schweitzer, Physica A191, (1992) 406. 

64. M. H. Cohen, E. N. Economou, C. M. Soukoulis, Phys Rev. Lett. 51, (1983) 1202. 

65. M. Schreiber, Phys. Rev. B31, (1985) 6146. 



